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Abstract 

The influence of the energy dependence of the free NN t-matrix on the op- 
tical potential of nucleon-nucleus elastic scattering is investigated within the 
context of a full-folding model based on the impulse approximation. The treat- 
ment of the pole structure of the NN t-matrix, which has to be taken into 
account when integrating to negative energies is described in detail. We cal- 
culate proton-nucleus elastic scattering observables for ^^O, ^''Ca, and ^''^Pb 
between 65 and 200 MeV laboratory energy and study the effect of the energy 
dependence of the NN t-matrix. We compare this result with experiment and 
with calculations where the center-of-mass energy of the NN t-matrix is fixed 
at half the projectile energy. It is found that around 200 MeV the fixed energy 
approximation is a very good representation of the full calculation, however 
deviations occur when going to lower energies (65 MeV). 

PACS: 25.40.Cm, 24.10-Ht 
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I. INTRODUCTION 



The scattering of protons and neutrons from nuclei has a long history as a tool 
to investigate the details of the reaction mechanism between a nucleon and a many nucleon 
system. The spectator expansion of multiple scattering theory is our theoretical ap- 
proach to define an optical potential for elastic nucleon-nucleus scattering. This expansion is 
predicated upon the idea that two-body interactions between projectile and target nucleons 
play the dominant role. 

In its most general form, the first order single scattering optical potential within 
the framework of the spectator expansion is given by the expectation value of the nucleon- 
nucleon (NN) transition amplitude and the ground state of the target nucleus. This 'full- 
folding' optical potential involves the convolution of the fully-off-shell two nucleon scattering 
amplitude with a realistic nuclear density matrix. In this form, the exact calculation of this 
optical potential requires a three-dimensional integration, in which the integration variable 
is coupled to the energy of propagation of the projectile and target nucleon. This very fact 
leads to a full-folding optical potential, which explicitly treats the off-shell behavior and 
the energy dependence of the NN t-matrix when carrying out the integration. Full-folding 
models along this line have been proposed and carried out with the conclusion that 
these model calculations based on the free NN t-matrix give a good description of the data 
for energies between 200 and 400 MeV, however are rather poor at 200 MeV and below. We 
have repeated these calculations within our computational framework based on a full-folding 
model employing a realistic nuclear density matrix and a free NN t-matrix where the off-shell 
behavior as well as the energy dependence are taken into account. However, we do not find 
as large effects as shown in Refs. 0,^. 

It can be argued that at intermediate energies the scattering of a projectile from a 
nucleon in the target nucleus may resemble free NN scattering mainly in forward direction. 
This is the justification for a common approximation to the optical potential in which the 
energy of the NN t-matrix is uncoupled from the integration variable by fixing it at half 
the projectile laboratory energy. We will study the accuracy of this assumption at different 
projectile energies. 

The structure of this paper is as follows. First we will review in Section II the 
relevant expressions for the single scattering optical potential in the impulse approximation 
as well as the full-folding procedure as used in our calculations. We will describe in some 
detail the numerical implementation for treating the energy dependence of the NN t-matrix, 
especially since this was left out in Ref. We then present in Section III elastic scattering 
results for proton scattering from a variety of nuclei in the energy regime between 65 and 
200 MeV, and end with concluding remarks in Section IV. 

II. THEORETICAL FORMULATION 

The transition amplitude for elastic scattering of a projectile from a target nucleus 
is given as 

Tei = PUP + PUGoiE)Tei, (2.1) 



2 



where P is the projector on the ground state \^a) of the target, P = ■'^^j^^, and Go{E) = 
{E — Hq + ie)~^. For the scattering of a single particle projectile from an A-particle target 
nucleus the free Hamiltonian is given by Hq = ho + Ha, where Ha stands for the target 
Hamiltonian. In the spirit of the spectator expansion the target Hamiltonian is viewed as 
Ha = hi + J2jf!^i Vij + H\ where hi is the kinetic energy operator for the ith target nucleon, 
Vij the interaction between target nucleon i and the other target nucleons j, and if* is 
an (A-l)-body operator containing all higher order effects. In a mean field approximation 
J2jjLiVij ~ Wi, where Wi is assumed to depend only on the zth particle coordinate. Thus, 
the propagator consistent with the first order in the spectator expansions is given as 

Go{E) ^ G,{E) = [{E - K) - ho -K~W, + ie]-\ (2.2) 

Here if*, having no explicit dependence on the ith particle, is replaced by an average energy 
f^*, which is at most equivalent to the separation energy of a nucleon from the nucleus. In 
most of our calculations we set i?* = 0, since the projectile energies are in comparison much 
larger. We will test this assumption. Since we are considering the optical potential in the 
impulse approximation, we neglect the mean field Wi of Eq. ( |2.2| ) in the following. 

The driving term of Eq. (|2.1j ) denotes the optical potential, which in first order is 

given as 

(k'|($A|P/7F|$^)|k) = t/(k',k) = (k'|($A|To.(^)|$A)|k). (2.3) 

i=n,p 

Here k' and k are the external momenta of the system, foi{£) represents the NN transition 
operator 

foi{E) = voi + Voigi{E)foi, (2.4) 

with 

g,{E) = [{E - E') -ho-h + ie]-\ (2.5) 



and Voi representing the NN interaction. The sum over i in Eq. ( |2.3| ) indicates the two 
different cases, namely when the target nucleon is one of Z protons, and when it is one of 
N neutrons. The energy £ is the energy of the interacting system. Inserting a complete 
set of momenta for the struck target nucleon before and after the collision and evaluating 
the momentum conserving 5-functions gives as final expression for the full-folding optical 
potential 

C/(q,K)= /rf'P^(P,q,K)fo,(q,^(^^K-P),£:) 

i=n,p 

Here the arguments of the NN amplitude foi are q = k' — k = /C' — /C and |(/C' + /C) = 
l(^K-P), where /C' = i(k'-(P-|-f ) and /C = i(k-(P + f-f ) are the nonrelativistic 
final and initial nuclear momentum in the zero momentum frame of the NN system, and 
K = i(k' + k). The factor ?7(P, q, K) is the M0ller factor for the frame transformation [§, 
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and Pi represents the density matrix of the target. Evaluating the propagator gi{E) of 



Eq. (2^) in the nucleon-nucleus (NA) center of mass frame yields for the energy argument 



E of the NN amplitude tqi of Eq. 



8 = E^j,- ^ ^ ^ • (2.7) 



4m 



N 



Here E]s[a is the total energy in the NA center of mass frame and mjy is the nucleon mass. At 
this point we assume ii^* = as was discussed earlier. Since we employ relativistic definitions 



for the NA kinematics we have Ej^ja = y (kc)^ + {mj^jc^Y + V (— kc)^ + (m^c^)^, where 
denotes the mass of the nucleus. With these definitions, the optical potential of Eq. 
becomes 



i=n,p 



r 1 yl + 1 i^^K. + Pr 

[/(q, K)=Y. J d'PviP, q, K) fo.(q, ^{^^ " P), E^A - ^ ^J^^ ) 



This expression shows that the evaluation of the full-folding integral requires the NN t- 
matrix not only off-shell but also at energies i?Ar^ > £^ > — oo. Specifically when going to 
negative energies, we have to take into account the pole structure of the NN t-matrix. The 
NN interaction supports a bound state at E^ = —2.225 MeV in the ^Si — ^ Di channel, the 
deuteron, and a virtual state in the ^Sq channel at -66 keV, the 'di-proton'. A virtual state 
means a pole of the NN t-matrix in the second sheet of the complex energy plane, which 
manifests itself on the real axis as a very narrow finite peak (about 100 keV wide) around 
the pole position. The deuteron, a true bound state, causes a pole in the NN t-matrix at 
the binding energy Ed- In order to explicitly treat this pole, we factorize foi into a pole term 
and a residue function 

To, = fo, X {E,,A - - ^^)' (2-9) 

4:171 N 



so that the optical potential takes the form 



U{q, K) = E / d'P 



i=n,p 



i,^K-P),i?.-.-(^l!^±^ 

iK+] 
4mjv 



Ena - ^ -Ed + te 



V{P, q, K)p,(P - P + )• (2-10) 

To obtain a simple pole, which can be treated with standard numerical methods, we perform 
a change of the integration variable to Q = ^^j^K + P, and Eq. ( p.lO|) becomes 



£>(q,K) = E J^Qf iQQ'- 4-. ' '^'■^"^ " ^' 



i=n,p 



r^(Q, q, K)p,(Q - ^(K + |), Q - - |)) 
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frin-D r ,nn^ a %(q> |(2K - Q), i?^^ - 
.f- J ^ Jo n2 _ n2 



i=n,p 



Qi-Q' 

viQ, q, K)p,(Q - ^(K + |), Q - - |)) 



-tYl JdQ 27rmNQd %(q, ^(2K - QdQ),Ea) 



i=n,p 



viQdQ, q, K)p,(Q,Q - + |), QdQ - - |)). (2.11) 



Here / represents the angular integration, Qd = \J Amiq{E]^A — Ed), and V denotes a 
Cauchy principal value integral. 



The cut in the integrand of Eq. ( 2.11 ) indicates the opening of a new channel in 



the scattering reaction. In this specific case it is the deuteron pickup channel, describing 
when an incoming projectile proton picks up a neutron and a deuteron is knocked out, thus 
removing flux from the elastic channel. In order to obtain some a priori estimate on the 
possible size of this additional channel on the elastic scattering reaction, it is worthwhile to 
look at experimental information on the deuteron pickup reaction. 

In the 1980's cross sections of the (p,d) reaction have been measured for the closed 
shell nuclei considered here [^|ll|] • For ^°Ca at 200 MeV, the differential cross section for the 
(p,d) reaction for small angles (10-20 degrees) is about 1 mb/sr while at 65 MeV it is about 
10 mb/sr. Compared to the typical size of differential cross sections for elastic scattering in 
a similar angular range (cp. Figs. 2 and 6), this is quite small. The effect of this additional 
channel will be studied in our calculations of elastic scattering observables and compared to 
calculations where the integration variable is decoupled from the energy of propagation of 
the projectile and target target nucleon. The latter is an approximation where the energy 
of the NN t-matrix is fixed at half the beam energy in the laboratory frame 

S = E, = -^ = -^-^^, (2.12) 
2 2m^ 2 2mN 

where kiab and are the on-shell momenta in the laboratory and NA system respectively. 
With this fixed energy approximation the full-folding optical potential of Eq. (|2.10|) becomes 



[/(q,K)= /rf'^^(P,q,K)fo,(q, 



i=n,p 



i(^K-P),.„ 



— iq^^ — iqN 

Full-folding calculations based on this approximation have been carried out by several groups 
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Fixing the energy at Eq as given in Eq. (|2.12| ) is a 'historic' choice, made by all 
calculations based on the KMT formulation [|13|. The argument for this specific choice of Eq 
is, that this is the energy at which free NN scattering at a fixed target nucleon would occur. 
This choice of Eq favors the forward scattering process, and it was argued that due to the 
narrow peaking of the density in momentum space, the scattering is dominated by forward 
scattering. Having in mind that the differential cross section for e.g. proton scattering 
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from ^"^Pb at 200 MeV falls off by five orders of magnitude between 5 and 30 degrees, this 
argument may capture some truth, however needs to be tested numerically. 

Before discussing the results of our calculations, we would like to elaborate some 
more on the details of our numerical implementation. Treating the pole structure of the NN 
t matrix while integrating over the energy of the propagation is standard in modern three- 



nucleon scattering calculations |Tj]. However, these calculations are carried out in partial 
waves, which allows treating the ^Si—^Di and ^5*0 channel separately. Our calculation of the 
full-folding optical potential of Eq. (|2.11| ) is carried out directly in three dimensions based 



on Monte Carlo integration. Thus, for the neutron-proton part of the optical potential, we 
have to treat the pole singularity and the virtual state simultaneously. The principle value 
integral of Eq. (|2.11|) is treated with standard subtraction techniques. When calculating 
the optical potential as given in Eq. ( |2.11| ), we have to interpolate the NN t-matrix in 
four dimensions, |q|, |2K — Q|, (2K — Q) • q, and the energy. For the momenta we use a 
three dimensional B-spline, for the energy a linear interpolation, since over a large range of 
energies the NN t-matrix is a slowly varying function of the energy. When carrying out the 
energy integration of Eq. ( |2.11| ), we find that for the higher projectile energies integrating 
out to -100 MeV cm. energy in the NN t-matrix is sufficient, while for the lower energies 
(65 MeV) we need to integrate out to -400 MeV cm. energy. 

Of course, the argument of the NN t-matrix being a slowly varying function of the 
energy is not true in the immediate vicinity of the virtual state in the ^5*0 channel. Since 
the peak around the pole position in the second energy sheet is finite and very narrow, we 
evaluate the integral over the t-matrix in this region separately and define 

1 /•Q0+A2 2K — Q Q'^ 

fav,i{ci, K, Ena) = -T — —T- / dQfoi{q, , Ena - ), (2.14) 

Ai + A2 JQo-Ai 2 ArriN 

where Qo is defined via E^a — 4;;^ = 0, which is close to the pole position on the second 
sheet of the complex energy plane. The 'average' fav,i is usually obtained from 40 energy 
points over a momentum interval of 0.5 MeV, and is then used as one of the interpolation 
points in the energy interpolation. The midpoint Qo for the integration in Eq. ( p. 141) depends 
via Ema on the projectile energy, which implies that the averaging process has to be carried 
out for each projectile energy. If we want to simplify the numerical procedure and create 
fav,i only once, we change the integration variable from Q to E and calculate 

1 pEo+5i 2K - O 
f,„,,(q, K, Eo) = — dEfo^iq, E), (2.15) 

ZOi JEo-Si 2 

where we have chosen £"0 = MeV and 61 = 0.25 MeV. The error in obtaining Tav,i from 
Eq. (p7[5|) instead of Eq. Eq. ( CT ) is small if is small. For our worst case tested, 



Ena = 50 MeV for proton scattering from ^osp]^ ^^le numerical error was 1%. 



III. RESULTS AND DISCUSSION 

In this paper the study of elastic scattering of protons from spin- zero target nuclei 
at energies that range from 65 to 200 MeV incident projectile energy is strictly first order in 
the spectator expansion and based on the impulse approximation. The full-folding optical 
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potential is calculated as outlined in the previous section, specifically as given in Eq. ( p. 111) . 
As a model for the density matrix for the target nucleus we employ a Dirac-Hartree calcula- 
tion |T^. The Fourier transform of the vector density, p{r',r), serves as our non-relativistic 
single particle density Another crucial ingredient in the calculation of the optical poten- 
tial, ?7(q, K), is the fully off-shell NN t-matrix. The calculations presented here employ the 
NN t-matrix based on the charge-dependent Bonn potential [0 . This potential is fitted to 
describe the Nijmegen database with a per datum ~1. It is also to be understood that we 
perform all spin summations in obtaining f/(q, K). This reduces the required NN t-matrix 
elements to a spin independent component (corresponding to the Wolfenstein amplitude A) 
and a spin-orbit component (corresponding to the Wolfenstein amplitude C). Since we are 
assuming that we have spin saturated nuclei, the components of the NN t-matrix depending 
on the spin of the struck nucleon vanish. The Coulomb interaction between the projectile 



and the target is included using the exact formulation as described in Ref. IT^ • 

At first we want to concentrate on proton scattering from different target nuclei at 
an intermediate energy. In Fig. 1 we display the differential cross section the analyzing 
power Ay, and the spin rotation function Q, for elastic proton scattering from ^^O. The solid 
line represents the full calculation of the optical potential according to Eq. (|2.11| ) and the 
dashed line shows the calculation where the energy of the NN t-matrix is fixed at half the 
projectile energy. It is remarkable how close the fixed energy result is to the full calculation. 
This may stem from a relatively weak energy dependence of the NN t-matrix and the fact 
that the folding with the sharply peaked density matrix does not sample the NN t-matrix 
in the negative energy region. We confirmed this numerically by artificially limiting the 
energy integration to positive energies and did not find noticeable effects in the scattering 
observables at 200 MeV. In Fig. 2 we show the observables for elastic proton scattering 
from ^°Ca at 200 MeV and in Fig. 3 those for elastic proton scattering from ^°^Pb. As 
in Fig. 1, the solid line represents the full calculation of the optical potential taking the 
energy dependence of the NN t-matrix into account, whereas the dashed line shows the 
calculation at fixed energy Eq given in Eq. ( p.l3| ). For all three nuclei the fixed energy result 
is remarkably close to the full calculation, a conclusion which essentially was also drawn in 
Ref. [0. A general trend is that for the fixed energy calculation the dip structure in the 
differential cross section and the spin observables is slightly more pronounced. 

In order to assess the importance of the additional energy shift given by E'^ in 
Eq. ( p.2|) , we show in Fig. 1 a calculation (including the energy dependence of the NN t- 
matrix), which was performed setting E^ = —8 MeV. This results in a shift of the total 



energy of the NA system to a slightly higher value, as suggested already in Refs. p8| , pJ 
The effect of this energy shift is negligible at 200 MeV as shown in Fig. 1. 

Next we turn to lower projectile energies, where we may expect to see differences 
between the full calculation including the energy dependence of the NN t-matrix and a 
calculation with a NN t-matrix energy fixed at half the projectile energy. In Fig. 4 we 
display the results for proton scattering from ^osp]^ -^qq ]\/[gY pjg 5 those for proton 

scattering from ^^O at 135 MeV. Compared to the results at 200 MeV, the difference between 
the solid line, which represents the full calculation, and the dashed line, which represents 
the calculation for a fixed energy of the NN t-matrix, becomes more pronounced, especially 
for Ay at larger angles as shown in Fig. 4. A similar figure is shown in Ref. [^], however a 
direct comparison between calculations is not possible, since different NN t-matrices as well 
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as densities are employed. The scattering observables given in Fig. 5 for proton scattering 
from ^^O at 135 MeV exhibit less difference between the solid and dashed lines, leading to 
the conclusion that even at this relatively low energy the fixed energy prescription in the 
optical potential is still amazingly and perhaps unexpectedly good. One common trend is 
becoming apparent in Figs. 4 and 5, namely that the differential cross sections predicted 
by the full calculations are systematically higher compared to the ones predicted by the 
calculations employing a NN t-matrix at a fixed energy. In Fig. 5 we also include a full 
calculation, in which the energy shift i?* is chosen to be = —8 MeV (dash-dotted line). 
Again, we conclude that the effect of this energy shift is minute. 

Last we want to consider scattering observables at an energy below 100 MeV pro- 
jectile energy. It has already been stated in the literature [|3|j5|j20(| that at these low energies 
the strict impulse approximation is insufficient to describe the experimental observables. We 
arrive again at the same conclusion, when showing in Fig. 6 the scattering observables for 
proton scattering from ^''Ca at 65 MeV. However, this is not the main point we want to 
make. When comparing a calculation including the energy dependence of the NN t-matrix 
in the optical potential (solid line) with a calculation where the energy is fixed at half the 
projectile energy (dashed line), we clearly see that the additional channels, the deuteron 
pick-up channel and the di-proton state have a sizable effect on the elastic observables. This 
is consistent with the experimental data, which show that the cross section for the (p,d) 
reaction is getting larger at lower energies. The differential cross section predicted by the 
full calculation (solid line) is higher than the one predicted by the calculation using the 
fixed energy prescription (dashed line), and actually slightly closer to the experiment. The 
spin observables obtained from the full calculations show much less structure compared to 
those from the fixed energy calculation. However, both impulse approximation calculations 
do not adequately describe the spin observables. A similar conclusion was drawn in Ref . |Q , 
although there the effect of including the deuteron and di-proton channel is more dramatic 
than in our calculations. For the first time we also observe a visible effect on the scattering 
observables, when we introduce an energy shift by setting = —8 MeV (dash-dotted line). 

From Fig. 6 we see that he impulse approximation, even if accurately and com- 
pletely calculated, is clearly inadequate in describing the elastic scattering observables at 
such a low energy, and additional effects have to be included. The first term in the specta- 
tor expansion contains as additional term the coupling of the struck target nucleon to the 
residual nucleus [If. For comparison we include a calculation in Fig. 6 as dotted line, where 
this term is calculated in an approximate fashion as described in Ref. [|r^. We see that 
this additional 'medium' contribution is necessary to get a better description of the spin 
observables. 

IV. CONCLUSION 

We have calculated the full-folding integral for the first-order optical potential using 
the impulse approximation within the framework of the spectator expansion of multiple scat- 
tering theory. The exact calculation of this full-folding integral requires a three-dimensional 
integration, in which the integration variable is coupled to the energy of propagation of 
the projectile and target nucleon. We have carried out the calculation taking into account 
the pole structure of the NN t-matrix when integrating to negative energies. Our optical 
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potentials are based on a Dirac-Hartree model for the nuclear density matrix and the charge- 
dependent Bonn potential for the NN t-matrix. Recoil and frame transformation factors are 
implemented in the calculation in their complete form. We calculate elastic scattering ob- 
servables for ^^0,^°Ca, and ^osp]^ projectile energies from 65 to 200 MeV laboratory 
energy and compare the full calculation with calculations in which the energy of the NN 
t-matrix is fixed at half the projectile energy. We find that this fixed energy prescription 
describes the full calculation remarkably well for proton scattering at 200 MeV projectile 
energy. This leads to the conclusion that the pole structure of the NN t-matrix does not 
play a role at intermediate energies. For projectile energies below 200 MeV we find that the 
influence of the deuteron and di-proton state slowly gains importance as we approach lower 
energies. However, between 100 and 200 MeV these effects are still relatively small. This 
is consistent with the small size of the experimentally measured cross sections for the (p,d) 
reaction. At 65 MeV we see considerable differences between the full calculation and the 
one in which the energy is kept fixed. In general, we find smaller differences between our 
full calculation and the fixed energy calculations as were suggested in Ref. 0. 

The calculations presented in this manuscript are strictly based on the impulse 
approximation using the free NN t-matrix. It has been stated in many places that the 
impulse approximation is insufficient at energies around 100 MeV projectile energy. We 
confirm this result, while carrying out complete calculations of the optical potential taking 
into account the energy dependence of the NN t-matrix. Additional corrections have to be 
taken into account at lower energies. Within the framework of the spectator expansion these 
corrections have been derived and carried out in an approximate fashion 0. A similarly 
complete calculation, as given here for the first order term in the impulse approximation, 
will require a complete Faddeev calculation of the correction term, which is at present not 
possible. 
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FIGURES 



FIG. 1. The angular distribution of the differential cross-section (^), analyzing power (Ay) 
and spin rotation function (Q) are shown for elastic proton scattering from ^^O at 200 MeV lab- 
oratory energy. The solid line represents the calculation performed with a first-order full-folding 



optical potential based on the DH density |15| and the CD-Bonn model [16| including the energy 
dependence of the NN t-matrix. The dashed line represents a calculation where the energy of the 
NN t-matrix is fixed at half the projectile energy. The dash-dotted line stands for a calculation in- 
cluding the energy dependence of the NN t-matrix and an additional energy shift by = —8 MeV. 
The data are taken from Ref. pl[ . 

FIG. 2. Same as Fig. 1, except that the target nucleus is ^°Ca, and the dash-dotted line is 
omitted. The data are taken from Ref. 



FIG. 3. Same as Fig. 2, except that the target nucleus is ^'^'^Pb. The data are taken from 



Ref. |2|. 



FIG. 4. Same as Fig. 3, except for ■^'^^Pb at 160 MeV proton kinetic energy. The data are taken 



from Ref. [g4|. 



FIG. 5. Same as Fig. 1, except for O at 135 MeV proton kinetic energy. The data are taken 



from Ref. [25 



FIG. 6. Same as Fig. 1, except for ^'^Ca at 65 MeV proton kinetic energy. The additional 
dotted line denotes a calculation including the effect of the coupling of the struck target nucleon 
to the residual nucleus. The data are taken from Ref. 12^. 
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